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Abstract 

Let  X{  be  i.i.d.,  EX{  ■■  0,  EXf  —  1,  Sn  •  X\  + ...  +  Xn,  f  be  a  reasonable 
function,  andTm  ■■  inf  |n  :  S»  >  m1/2#  (^)j-  The  limiting  distribution  of  $rm- 
m1/2#(^)  is  found,  and  used  to  compute  heuristic&lly  higher-order  correction 

terms  to  the  passage  probability,  P  j^p  <  tj,  the  first  term  of  which  is  given  by 
the  invariance  principle. 

Key  Words:  Random  Walk,  Diffusion  Approximation,  Passage  Times,  Tests  of 
Power  I. 


Much  of  the  motivation  and  interest  of  the  problem  considered  in  this  paper 
arises  where  the  underlying  object  of  interest  is  the  time  at  which  a  random 
walk  exits  from  a  certain  region.  Typically  this  cannot  be  calculated  exactly 
and  some  sort  of  approximation  is  necessary,  which  proceeds  by  imbedding  the 
boundary  of  the  region  in  an  infinite  family  of  similar  boundaries,  allowing  them 
to  move  out  to  oo  and  calculating  the  distribution  of  the  exit  time  asymptotically. 
Imbedding  the  boundary  in  a  family  of  boundaries,  by  the  method  employed  here, 
is  equivalent  to  imbedding  the  random  walk  in  a  family  of  random  walks,  and 
considering  the  time  at  which  they  cross  a  fixed  boundary.  This  can  be  done  in 
such  a  way  that  the  imbedded  random  walks  tend  to  become  degenerate  along 
the  line  of  drift,  and  thus  exit  the  region  in  a  small  neighborhood  of  where  the 
line  of  drift  exited.  This  means  that  only  the  properties  of  the  boundary  at  the 
point  at  which  the  line  of  drift  crosses  affect  the  distribution  of  the  excess  over  the 
boundary  and  the  crossing  time.  This  seems  to  be  a  very  bad  property,  because 
boundaries  which  are  entirely  different  but  happen  to  be  tangent  at  one  point 
appear  the  same  to  random  walks  which  tend  to  cross  at  that  point,  and  one 
fears  that  the  asymptotic  calculation  has  misled.  The  monograph  of  Woodroofe 
(1082]  discusses  this  case. 

An  alternative  technique,  to  which  the  above  objection  cannot  be  made,  is 
to  scale  the  random  walk  so  that  it  converges  to  a  Wiener  process,  possibly 
with  drift.  Then,  the  probability  that  the  random  walk  crosses  a  boundary  in 
a  given  interval  of  time  approaches  the  probability  of  the  same  event  for  the 
Wiener  process.  The  main  objection  that  one  can  make  to  this  calculation  is 
that  anything  that  distinguishes  one  distribution  from  another  is  wiped  out  in 
taking  the  limit.  The  proposed  remedy  is  to  calculate  higher  order  correction 
terms  to  the  Wiener  process  crossing  probabilities  which  are  typically  highly 
distribution  dependent. 

*n  jh'*  du 

Section  Intakes  up  the  question  of  finding  the  distribution  of  excess  over  the 
boundary  for  random  walks  normalized  to  converge  to  a  Wiener  process.  In  some 
cases,  such  as  tests  of  power  1  which  are  discussed  in  Section  2,  the  connection 
of  the  excess  over  the  boundary  to  the  problem  of  finding  the  correction  to  the 
Wiener  process  crossing  time  is  simple  and  direct,  and  in  other  cases,  especially  - 


that  of  Siegmund  [1979],  quantities  connected  with  the  excess  over  the  boundary 
appear  in  ways  that  do  not  shed  much  light  on  the  general  problem.  Section  3 
gives  a  heuristic  method  for  correcting  the  random  walk  probabilities,  based  on 
a  reflection  principle  method  of  Siegmund  and  Yuh,  which  is  the  author’s  best 
idea  for  the  connection  between  the  excess  over  the  boundary  and  the  corrected 
diffusion  approximations  to  the  hitting  times  in  the  general  case.  The  results  of 
Section  1  are  used  here. 


I.  Excess  Over  the  Boundary  in  Diffusion  Approximations 


Here  is  some  notation  that  will  be  used  in  Theorem  1.1.  Let  y+(t)  ,*-(*)  >0  < 
I  <  1  be  continuously  differentiable  functions  such  that  g+  >  0,  <  0.  Let  Sn 

be  a  nonarithmetic  random  walk  with  E|Si|3<  oo,  ESy  —  0,  and  ESf  l.  Let 
Tm  —  inf  {«  <  m  :  Sn  >  m1/2#*  (&))  or  Sn  <  to1/2#-  (&)}.  Let 

J*+(r«i)  —  («Tm  -  ™'/sl+  (~))  >(Srm  >0)  ’ 

and  be  a  random  variable  with  the  distribution  of  excess  over  a  distant 
boundary.  That  is,  if  ra  —  inf  (a  :  Sn  >  «},  then  Su  -  a  The  conver¬ 

gence  in  distribution  is  a  consequence  of  the  renewal  theorem.  Ta  denotes  the 
field  generated  by  {X,l(a>f)},  where  X,-  ■■  5,*  -  5,_j.  L(X)  will  denote  the  dis¬ 
tribution  of  the  random  variable  X,  and  (<1,(2)  C  [0,1].  Let  VT(<)  be  a  standard 
Wiener  process,  and 


T  -  inf  {t  <  1 :  W(t)  >  f+(t)orW(t)  <  #-(*)}. 


Theorem  1.1.  With  notation  as  above,  and  ((1,(2)  C  [0, 1], 

p{“  e  (ti,t2),sTm  >  o,f?+(rm)  >  *}  -  P{T  6  (*1,*2),  W(0  >  0>P{*(Si)  >  r). 

Furthermore,  if  E\S\\^  <  oo,  1  >  3  then  there  is  a  random  variable  Z  >  0  with 
EZ7-3  <  00i  and  P{Sfm  >  0,E+(rm)  >  y,  Tm/m  <  1}  <  const.P{£  >  y). 

The  facts  about  the  hitting  time  itself  are  a  standard,  though  not  completely 
trivial,  application  of  the  invariance  principle.  Specifically 


—inf  (n  :  Sn  >  m1/2  g  (— )}  -■  inf  (— 
m  l  w)  »  im 


<  1 


m 


172 


>f 


0} 


—  inf {t  <  1 :  W(t)  >  f(«)>. 


S 


It  is  therefore  taken  as  given  that  -*  T.  The  key  to  establishing  Theorem 
1.1  is  Lemma  1.1,  which  implies  that  even  as  the  boundary  moves  out  to  oo  and 
smooths  out,  the  process  can  get  within  k  units  of  it,  where  k  is  some  positive 
constant,  without  (with  high  probability)  crossing.  The  proof  is  then  completed 
by  observing  that  starting  from  reasonably  close  to  the  boundary  the  process 
crosses  as  if  it  were  crossing  a  distant  constant  boundary. 

Lemma  1.1.  The  0(1)  Boundary  Crossing  Lemma.  With  notation  as  in 
Theorem  1,  there  is  a  k  >  0  and  a  random  variable  Z  >  0  with  EZ 1,-3  <  oo 
depending  on  the  boundary  only  through  a  —  8aPtg|o,ij  f+(0>  such  that 

P(R+(rm)  >  <  1)  <  kP{Z  >  ,). 

m 


Proof.  The  distribution  of  excess  over  the  boundary  can  be  written 
PfR+ITm)  >  r,  5=  <  1}  -  £  p{sH1  >  m>/2  f+(j±i)  +  ,,  Tm  >  /} 

-  '■k.  -  *  >  ™,/i  m^))*  » + *} 

X  p{5j  €  m1/2 -  **,  Tm  >  ;] 

</  p{Sj+l-Sj  >!  +  »  +  *}  (1-1) 

X  €  m1/2  f(^)-dz,  Tm  >  >} 

<  ]£  plsHl  -Sj>n  +  v) 
n«l 

x  Jr, p  iSi  6  ( m'11  •*  (™) "  *•  m‘/2  •*  »}  - 

where  the  next  to  the  last  line  follows  from  m1/2|?+  (^)> 
for  m1/2  large.  It  is  necessary,  therefore,  to  bound 

*Vn-Y*P {5i  €  (m1/2 1+  (£)  -  n,  m1/2  f+  (^)“n  +  1)>  T™  >  /} • 

This  is  the  expected  number  of  visits  to  a  certain  set  by  the  space-time  random 
walk  before  exiting  from  a  larger  set.  Assume  that  the  random  walk  has  the 


4 


property  that  P{Sn *  >  (1  +  <*)n}  >  t  >  0  Vn.  Assume  that  n  <  m1/2.  Then 


F{V„  >  *  +  n2|V«  >  *}  <  p|sn,  <  a  -j^  +  nj 

<WB*<(1  +  ®W 

<  1-T 

The  first  inequality  merely  bounds  the  probability  of  spending  n2  more  time 
units  in  the  set  before  exiting  the  larger  set  by  the  probability  that  in  precisely 
n2  more  steps  the  random  walk  has  exited  the  larger  set.  It  is  easy  to  see  that 
this  implies  EVn  <  en2,  where  the  e  depends  only  on  7,  which  depends  on  the 
boundary  through  a  only.  Also,  because  7m  <  m,  EVn  <  en2  is  very  obvious 
for  n  >  m1/2.  To  remove  the  restriction  imposed  above,  note  that  the  CLT 
guarantees  that  the  inequality  holds  for  n  >  hq  and  the  entire  argument  can  be 
repeated  on  P{Vn  >  k  +  non2|Vn  >  k)  to  produce  a  bound  EVn  <  en2,  where  e 
could  depend  on  no-  So  the  bound  EVn  <  en2,  where  e  depends  on  the  boundary 
only  through  a  has  been  established.  Substituting  this  bound  into  (1.1)  gives 


p{r+(T„)  >  »,  ^  <  l}  <  «  Z  »2  />{«,>„  +  »). 

'  m  1  m-1 

Now  if  F|5j|3<  oo,  Z  can  be  chosen  to  be  the  random  variable  with  distribution 

Pl ,  ^  _  Em. l  *2P{*1  >  *  +  *») 

P{Z>,) 

If  i  <  3  the  statement  about  the  moments  follows  from 


oo 

EZ^z  <  const.  E  p~4  P{Z  >  j) 

o 

00  00 

—  const.  £  p~4  5;  n2  P{Si  >  j  +  n) 

J—0  n- y 

«-  const.  £  £/y_4n2P{5|  >/  +  n) 
n— iy— 1 
00 

<  const,  £  n7-1  P{5|  >  n) 

B"1 

<  00, 


5 


since  £|5iP<  oo,  and  where  the  constant  is  not  necessarily  the  same  from  line 
to  line. 


Proof  of  Theorem.  Let 

-inf  jn  <  m-Sft  >  m1/2 y+  - k or Sn  <  m1/2 f_ + *J , 

where  it  is  assumed  that  the  relation  beteen  k  and  m  is  such  that 
m1/2 #+  (%)-k>  0,  rn1/2 i- (&)  +  *  <  0.  Denote  by 

*•<**>-  (sr*  -  (m,/J  »+  (£)  -  *))  '(Sr„  >0). 

and  similarly  for  R~(Tj^).  Lemma  1  applies  to  R±(t£1),  so,  t  >  0  being  given,  it  is 
possible  to  choose  k  so  that  P  {R±(Tjil)  >§}<«.  It  may  also  be  assumed  that  k  is 
large.  Fix  an  interval  (tj,*2).  The  first  step  is  to  show  that  those  paths  for  which 
Tm  €  (mf  1,01*2)  and  STm  >  0  are  essentially  those  for  which  6  (mf1,mt2)  and 
STk  >  0.  The  symmetric  difference  between  the  two  sets  is  contained  in 

{T&  €  (mti,mt2),  Sj*  >  0,Tm  €  (m*1,mt2),  STm  <  0} 

U  {T*  €  (mli,m*i),  Tm  g  (mfi,ro*2)}  ^ 

U  {Tm  €  (m*i,m*2),  T*  g  (ml,,  m*2)} 

U  {T&  6  (m*i,m*2),  Sy*  <  0,  Tm  6  (m*i,m*2),  STm  >  0). 

Consider  the  first  set.  Let  <1  —  sup{y_(f),t  €  (<i> <2)}*  f  ■*  inf{y+(*M  €  (*i,<2))- 
By  hypothesis  u  <  0  <  f.  By  the  strong  Markov  property 

P{7£  €  (m*i,ro«2),  Sy *  >  0,  Tm  €  5Jfn  <  0} 

where  r/  —  inf{(  :  S*  (J  /},  and  P*  denotes  the  measure  governing  the  random 
wlk  started  from  the  point  z  at  time  0.  This  is  so  because  within  logm  after  Tm 
it  exits  when  the  upper  boundary  is  within 


of  the  value  of  the  random  walk  at  time  and  where  0  <  X  <  1.  It  is  clear 
that  the  probabilities  of  both  the  first  and  the  second  events  go  to  0  as  m  -*■  oo. 

The  second  set  in  (1.2)  b  also  easy  to  handle.  Similar  to  the  above 

P{Tm  6  (ml|,mf2),rm  £  (m(ltm<2)} 

<  P{r£  6  (mt2  -  m1/2,  mt2)> 

+  p*<r(-oo,t+i«)  >  ™I/2>  +  '*<*•(- >  •»,/2>- 
It  b  easy  to  show  that  the  exit  dbtribution  of  the  Wiener  process  through  bound¬ 
aries  considered  here  b  continuous.  Consequently,  the  invariance  principle  shows 
that  the  first  probability  goes  to  0,  and  the  second  and  third  approach  0  as  above. 
(Strassen  (1067]  shows  that  in  the  one  boundary  case,  the  exit  dbtribution  actu¬ 
ally  has  a  continuous  density.  Hb  argument  presumably  extends  to  the  two-sided 
case,  but  thb  strong  a  result  b  not  necessary  here.) 

The  last  two  probabilities  can  be  handled  analogously. 

Therefore,  the  above  argument  shows  that  regardless  of  k,  for  m  sufficiently 
large, 

P({Tm  €(m*i, mt2),  STm  >  0,  R+{Tm)  > 

A{T*  €  (ml|,mt2),  Sr*  >  0,i?+(rm)  >  *})  <  e/2, 

where  e  >  0  b  arbitrary.  So  for  k  sufficiently  large,  Lemma  1  guarantees  that 
P{{Tm  €(ml|,wi<2),  Sj>m  >  0,  R+(Tm)  > 

A (4.  nR+{Tm)  >  *))  <  c, 

where  A^  —  {T&  €  (m<i,mt2),  Sj.*,  >  0 ,P+(7&)  <  ^}.  The  next  part  of  the 
argument  will  show  that  conditionally  on  7Th  ,  on  A^,  exiting  the  region  b 
much  like  crossing  a  distant  constant  boundary.  Some  extra  notation  must,  alas, 
be  introduced  to  show  thb. 

If  Zn  b  any  random  walk  with  0  <  EZ\  <  oo  denote  by  ta  *■  iuf{n  :  Zn  >  a), 
Ra{Zi)  —  and  Roo{Z\)  a  random  variable  with  the  dbtibution  L(Roo(Zi))  — 

lima->oo  L(Ra(Z)),  which  limit  exbts  by  the  renewal  theorem.  The  starting  point 
b  the  following  identity: 

P(4.  A  <P+(rm)  >  r)  -  f  P{P+(rm)  >  »|/r*  })dP. 


(1.3) 


Consider 


P{R+(Tm)>  v\?Tk) 

“  P{R+{Tm)  >  v,Tm-Tm  <  logm|7r>*  }  (1.4) 

+  P{R+(Tm)  >  y,  Tm  —  T%l>  logm|^i  }. 

As  in  the  proof  of  the  bounds  in  (1.2),  the  second  probability  goes  to  0  for  m 
sufficiently  large,  regardless  of  4.  To  bound  the  first  probability,  note  that  over 
an  interval  of  length  log  m  the  boundary  varies  little,  that  is 


sup  m 
j  <  log  m 


logm 

m1/2 


0. 


Let  us  pick  an  and  compare  the  value  of  Tm  to 

-  W {>  >  Tt, :  Sj  - sr»  >  - sTj,  ± ,,} . 

On  the  set  A&,  w*1/2  §  -  Sj*  >  4/2.  Suppose  4  is  chosen  so  that 

•nf#>£/4 P{Ro(S\)  >  2<|}  >  1  - Then  on  A*, 

PiTi,  >  KVl *)  -  A*  {*m„>  <  2«,|7X.}  <  <■ 

So  since  >  Tm*1,  P{Tt £,  —  Tm1  |7y*  )  >  1  -  e  on  A^,.  Also  it  is  clear  that 

4  HI 

{7#  -  I'm*1.  7m  -  T*  <  logm}  C  (7#  -  rm), 


and  further,  that  the  excess  over  the  three  boundaries  differs  by  <  2<j  on  this 
set.  Therefore,  choosing  4  so  that 


sup  \P{Ra  >  y  +  (i)-P{Ra  >  *r  —  «i)I<  «> 
o^*/4 

for  m  large  the  first  term  on  the  right  hand  side  of  (1.4)  can  be  bounded  above 
and  below  as  follows: 

P{Wml)  >  y  +  2e,  7#  -  T*  <  logml^J  -  P{T%  *  Tm\hm) 

<  P{R(Tm)  >V,Tm-T*l  <  logrn|^m} 

<  P{R(T%)  >  y,  T%  ~Tm  <  \ogm\frJ  +  P{TZ(  +  Tm\^J. 


Considering  only  the  upper  bound,  it  in  turn  can  be  bounded  by 
PWTM  >  irl^J  +PlTt,riTm\rrj. 

For  m  large  enough,  the  second  term  can  be  made  <  $  on  A&,-  Also,  Jb  can  be 
chosen  sufficiently  large  so  that 


Then  on  a£, 


sup  \P{Ra(Si)  >  y)  -  P{R<x>{Si )  >  y}|<  c. 
a>k/4 


\P{R(T%)  >  y\*rj  -  P{Roo(Sx)  >  y}|<  e. 


Substituting  these  various  bounds  into  (1.4)  gives 

P{Tm  €  (mtltmt2),  STm  >  0,  R+{Tm)  >  v)  <  P{A^n(R(Tm)  >  y)}  +  e 

<P{A*l)P{R00>y}  +  it 

<  P{r£  €  (mt|,m«2),  ST *  >  0}P{Roo  >  y}  +  5< 

<  P{T  €  WW  >  Oj/MKoo  >  y)  +  <k 

for  m  large  by  the  invariance  principle.  The  lower  bound  can  be  established 
similarly,  completing  the  proof. 

Remarks. 

1)  The  one-boundary  case  follows  from  the  two-boundary  case  by  letting  the 
lower  boundary  (say)  -» — oo. 

2)  The  restriction  y+  >  0,  #_  <  0  is  not  necessary,  it  can  be  replaced  by 
y_(t)  <  y+(«),  0  <  t  <  1.  It  is,  of  course,  necessary  to  assume  that  y_(0)  < 
0  <  y+(0),  because  otherwise  the  Wiener  process  crosses  immediately  at  0  (or  the 
lower  boundary,  say,  is  actually  an  upper  boundary). 

3)  A  drift  can  be  added  to  the  random  walk  at  stage  m  to  make  it 
converge  to  a  Wiener  process  with  drift,  and  this  is  equivalent  to  adding  a  linear 
term  to  the  boundary.  In  particular,  the  excess  over  the  boundary  is  the  same. 
Other  ways  of  getting  a  Wiener  process  with  drift  are  possible,  e.g.  at  stage  m 
choose  the  member  of  a  fixed,  one  parameter  exponential  family  with  expectation 


—^2  to  generate  the  random  walk.  The  main  new  tool  needed  for  this  case  is  a 
uniform  rate  of  convergence  for  Ra  to  Roo,  over  members  of  the  exponential  family 
in  a  neighborhood  of  the  distribution  with  drift.  In  an  exponential  family,  it  is 
not  hard  to  show  (possibly  with  the  additional  assumption  of  nonarithmeticity) 
that  i?oo(Sf)  is  a  continuous  function  of  I.  Therefore,  if  for  large  a  the  quantities 
Ra  for  the  different  distributions  are  all  close  to  their  limiting  value,  they  must 
all  be  close  to  the  value  of  Roo  for  the  mean  zero  member.  It  is  believed  that 
an  argument  of  this  sort  can  replace  the  part  of  the  proof  following  (1.4),  and 
that  this  can  be  made  rigorous  using  a  coupling  argument,  but  details  have  not 
been  worked  out.  Siegmund  and  Yuh  [1082]  consider  this  framework.  They  use 
a  different  method  of  establishing  this  uniform  convergence. 

4)  There  is  a  clear  conjecture  for  what  happens  in  the  lattice  case.  Con¬ 
sider  the  case  of  a  single  upper  boundary,  denoted  by  g.  Suppose  that  Sl  is 
arithmetic  with  span  1.  The  technique  of  the  proof  above  shows  that  the  dis* 
tribution  of  R+{Tm)  is  approximately  floo(Si)  -  (m1/2  g  (7^)),  where  (*)  =  *- 
[x]  =  fractional  part  of  x,  and  where  -Roo(Sj)  is  independent  of  Tm.  The  random 
variable  g  [  if)  is  converging  in  distribution,  and,  if  g  is  not  constant,  the  lim¬ 
iting  random  variable,  g(T),  has  a  reasonably  smooth  density.  When  X  is  any 
bounded  random  variable  with  a  continous  density  it  is  not  hard  to  show  that 
(m1/2  X)  U,  where  U  is  a  uniform  random  variable  on  [0, 1],  Therefore  the 
obvious  conjecture  is  that  f?+(rm)  -*  Roo{S\)—U,  Roo{S\),  U  independent.  There 
is  a  problem  with  rates  of  convergence,  or  perhaps  local  uniformity  in  time  of 
convergence  that  prevents  this  argument  from  being  rigorous. 

2:  Application  to  a  Test  of  Power  I. 

The  reader  shoule  recall  that  the  main  motivation  for  this  work  is  the  cor¬ 
rection  of  the  power  functions  of  statistical  tests  defined  by  boundary  crossing. 
It  is  assumed  that  the  crossing  distribution  for  the  Wiener  process  is  known,  or 
very  well  approximable  and  the  problem  is  to  correct  this  crossing  probability  for 
discreteness  or  non-normality.  The  following  is  an  example  where  information 
about  the  excess  over  the  boundary  alone  suffices  to  obtain  the  first  correction 
term  of  the  Wiener  process  boundary  crossing  probabilities.  The  crossing  prob¬ 
lem  below  is  corrected  with  tests  known  as  tests  of  power  1,  and  are  discussed  in 
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Woodroofe  [1082],  Chapter  6. 


Let  Pfi  be  a  family  of  probabiity  measures  under  which  the  process  {St;  t  > 
0}  is  a  Wiener  process  with  mean  jit  and  variance  t.  Define  a  measure  Q  by 

Q(A)-  f 

where  Q  is  a  standard,  normal  distribution.  At  time  t,  Q  <<  PQ  and  the  likeli¬ 
hood  ratio  is  given  by 


Lt 


*Q_ 

dP0 


i1_(riT),/Vw<i«>. 


Let  e  >  1  and  define  a  stopping  time  T  by 


T  -  inf  {t :  Lt  >  e) 
«=inf{t:|Sf|>  f(t)) 

where 


f(0  —  (2(t  +  lXloge  +  log(l  + 1)))1/2 . 

The  fundamental  identity  of  sequential  analysis  (a.k.a.  Wald’s  likelihood  ratio 
identity,  see  Siegmund  [6]  propositions  2.24  and  3.2)  shows  that 

Pq{T  <  oo}  *=>  T  <  oo^ 

1 

e  ’ 

since  P^T  <  oo}  —  1  V#i  yd  0.  The  discrete  version  of  this  test  defines  a  stopping 
time 


Tm  “  inf{n  G  Z+;  Lnjm  >  e}f 

which,  because  of  the  scaling  properties  of  the  Wiener  process  can  be  written  as 
being  equal  in  law  to 


Theorem  3.1. 


P0{Tm  <  oe}  -  i  -  -^^(^(S,)) 

xM*|Ti)i:r 


T  <  oo 


l-U) 


Proof.  By  the  fundamental  identity 


Po{Tm  <  oo}  —  Eg  ; 

-E*kk) 


Tm  <  00 


-  Ill B-  {"»  -  (nffffil)**  +  itwI-Ml)}  «' 

+'(^) 


4>{»W 


where  R±  is  the  excess  above  or  below  the  boundary.  In  this  scale  one  expects 
R±  —  0  Each  of  the  summands  in  the  brackets  in  (2.2)  is  non-negative, 

and  R±  is  always  a  positive  quantity.  A  one  term  Taylor  expansion  gives  that 


P0{Tm  <  oo}  " 


tjTmfm) 
2(7m/n*  + 


where  fm(  )  is  &  random  quantity  satisfying  |/m(-)l<  1  and  fm{  )  -*  1  as  m  -►  oo. 
It  is  clear  how  the  integrand  should  behave.  Theorem  1.1  gives  the  asymptotic 
distribution  of  m1/2  R±,  and  the  fact  that  it  is  asymptotically  independent  of 
is  and  continuous  bounded,  so  the  first  summand  is  converging  in  distribution, 
and  the  second  is  converging  to  0.  The  only  fiy  in  the  ointment  is  that  Theorem 
1. 1  does  not  give  quite  enough  information  to  claim  uniform  integrability  for 
| m1/2  /?±|  or  J m1/2  A±|2,  which  also  must  be  uniform  over  p.  This  must  be  taken 
care  of  by  a  special  argument. 


Under  P ^  Tm  can  be  written 

-w{";|Sn"^T72' n|>  m‘/2«D}  ’ 

where  Sn  is  a  N{0, 1)  random  walk.  That  is,  Tm  is  the  sort  of  stopping  time 
considered  in  Theorem  1,  with  drift  Similarly,  R±  can  be  written  as  being 
equal  in  law  to 


The  technique  is  very  much  like  the  proof  of  Theorem  1.1.  Suppose  one 
wants  to  establish  bounds  on  the  distribution  function  for  excess  above  the  upper 
boundary.  This  distribution  can  be  written 

P.{Rm>v)-ftP{S1>i  +  z) 
and  consequently 

P{Rm>j)<  £  P{Sx>i  +  n) 
m— 0 

X  £p.{s,€m,/2,(^)-<-,n  +  l)}. 

as  in  Theorem  1.1.  If  the  bound  of  const,  n5  is  obtained  on  the  second  sum, 
the  Rm  has  been  stochastically  bounded  by  a  distribution  with  moments  of  all 
orders,  as  in  Theorem  1.1.  Therefore,  we  seek  to  control 

B  £  ‘o  j+i)  (s" + ”,I/a » 
in  the  case  of  the  lower  boundary  or 

E  £  iys+i)  (m!/2 1  (^!r)  “  5") 


I 


(2.4) 


in  the  case  of  the  upper.  There  are  naturally  two  cases,  the  drift  0  case  which 
works  equally  well  for  small  drifts,  with  |/fc|<  *o,  and  the  case  of  large 
drifts.  In  the  zero  drift  case  the  upper  and  lower  boundary  are  the  same,  but 
in  the  positive  drift  case  this  is  not  so.  It  will,  however,  be  clear  in  the  positive 
drift  case  that  the  proof  of  the  zero  drift  case  works  equally  well  with  respect  to 
boundedness  above  the  upper  boundary,  so  only  excess  below  the  lower  boundary 
need  be  considered.  The  case  of  a  sizeable  negative  mean  follows  by  symmetry. 

Case  1.  Zero  Drift.  The  argument  for  Theorem  1.1  provides  a  bound  on 
(2.3)  up  to  some  time  (gm  and  also  for  visits  within  k^mM2  of  the  boundary. 
Therefore  assume  j  >  k^m1!2.  Note  that 

i,  .  loge  +  l+log(l  +  *)  1  /log(l  + 

“  2((1  +  zXlof  *  +  Ml  +  s)))1/2  ~  2V  *  ) 

Suppose  that  n  >  ;'5.  Then 


P(exiting  within  j5  more  stepsjwithin  j  of  boundary  at  time  ti) 


(2.5) 


But  it  is  easy  to  see  that 


and  hence  by  scale  properties  of  the  normal  distribution  the  quantity  in  (2.4) 
is  bounded  away  from  0.  This  supplies  a  bound  of  const.  (/*)  on  (2.3),  and  the 
proof  of  the  claim  is  finish  i  a-  in  Theorem  1.1.  Note,  as  explained  above,  that  if 
the  random  walk  had  beu  a  %  upwards  this  would  have  bounded  the  excess 
uniformly  above  provided  >h  i-.e  mean  drift  is  bounded  above,  which  is  only 
to  say  the  increments  are  sto  hastically  bounded. 


Case  2.  Positive  Drift.  Let  p  —  ^g,  k  >  1.  Consider  j  in  two  classes, 


>  <  }  and  j  >  J.  For  j  >  J 


I 


v;,f 


m 


■:1 


yf 

:4 


:a 


si 


ss 


'4 


1 


i 


£ 


fi.{s„+jj  >  -m'/2*(“T“)  +  y+  1|S»  €  »(1^j + U.i + 1)1 
>  PpiSp-j2!*  >  2j+  1}  >  t  >  o  vy. 


And  by  Lemma  2  below 


P{ever  being  within  j  of  boundary  [farther  than  2 j  from  boundary)  <  t  . 


This  provides  a  bound  on  (2.3)  in  the  manner  of  Theorem  1.1.  On  the  other 
hand,  for  j  <  —  i,  since  is  bounded 


<  -m*/s  ,  )  ]S„  +  m1/2  ,  (£)  +  (/,/+ 1)) 


But  ^  <  1,  the  left  hand  side  of  the  inequality  is  a  N[0,j2)  and  the  right 

hand  side  is  (uniformly)  within  its  standard  deviation.  This  gives  a  bound  on 
(2.3)  of  const,  (y2)  as  above,  and  finishes  the  proof  of  the  claims. 

Now  it  has  been  shown  that  the  quantities  represented  by  m}l2R±  in  (2.3)  are 
uniformly  integrable  for  all  values  of  the  parameter.  To  finish  the  proof  of  the 
theorem  it  is  a  simple  matter  to  take  the  limit  in  (2.3)  and  using  the  asymptotic 
independence  of  mxf2R±  and  Tm,  and  the  form  of  the  limiting  distribution  given 
in  Theorem  1.1,  to  write 


/wrm  <  =o)  -  i{i + /  ”  E>  (_y } 


where  the  definitions  of  Q  and  of  the  stopping  time  have  been  used.  This  finishes 
the  proof  of  the  claim. 


Lemma  S.l.  Let  p  >  0,  z  >  0,  zq  >  0,  and  rx  —  inf{n  :  Sn  >  a}-  Then 

P-p{rxt/it  <  °°}  <  «"2*0- 


Proof. 

P-pirx0/p  <  °°}  “  ^/i{exp(-2/i(Sr>o/(|)} 

”d  s%/. 

The  form  of  the  approximation  in  Theorem  2.1  illustrates  the  main  theoret¬ 
ical  and  practical  drawback:  it  depends  globally  on  the  Wiener  process  cross¬ 
ing  distribution.  Practically  this  means  that  to  use  this  method  it  would  be 
necessary  to  compute  the  crossing  distribution  for  the  Wiener  process,  which 
is  usually  impossible  to  write  down  explicitly.  Theoretically,  related  questions 
such  as  P0  <  1 J  necessitate  having  to  approximate  Pg  <  l  j  which  does 
not  involve  simply  the  excess  over  the  boundary.  A  method  for  approximating 
Pofy  <  1  j  f°r  certain  random  walks  and  boundaries  is  the  subject  of  Section 
3* 

3.  Corrected  Diffusion  Approximations  to  First  Passsge  Times. 

First,  some  notation  will  be  introduced  and  summarised.  It  will  be  assumed 
that  the  probability  space  will  accomodate  a  random  walk,  Sn,  with  ES\  —  0, 
ES$  —  l,  ES%  —  0,  and  a  standard  Wiener  process  W(i).  T’s  will  be  used  for 
random  walk  crossing  times  and  r’s  for  those  of  the  Wiener  process.  rg  will  denote 
different  thinp  in  the  one-sided  and  two-sided  case;  the  reader  may  consider 
himself  fairly  warned.  It  will  generally  be  clear  what  boundary  the  random  walk 
is  crossing,  and  Tm  will  denote  that  crossing  time  for  -%j.  If  t0  —  inf  {n :  Sn  > 
•},  then  it  is  known  from  renewal  theory  that  -  a  converges  in  distribution 
and  is  uniformly  integrable.  Set  $  —  lima-*oo  E  ST,-a.  The  results  are  as  follows: 

(I)  Let  f  <  0  be  continuously  differentiable  on  (0, 1],  t  <  1, 

Tm  —  inf  |n  :  Sn>  m */2  f  (^)J» 
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and  tj  mm  inf{t :  W(t)  >  f(t))  for  any  function  /.  Then 


and 

w*  <  '•  "W  <  •«>  -  *)  +  • 

(2)  Assume  that  X\  is  symmetric.  Let  §  >  0  be  continuously  differentiable  on 
[0,1],  0  <  <0  <  *1  <  1,  Tm  -  inf{n  >  mt0  :  |Stt|>  m1/2*^)},  and  rh  -  inf{t  > 
<o :  |Wr(r)|  >  h(t)}  for  any  function  h.  Then 

+  P{t°  K  Tg+0/m'l*  <  *1>  + 
and 

P{<0  <  ^  <  *!’  *"1^2(»(ll)-*)} 

—  /'{to  <  Tg+p/mi/*  <  *l>  |W(*l)l<  »(*l)  -*}-♦-  o  (^172)  * 

Notice  that  because  of  the  assumed  symmetry  of  the  random  walk,  to  consider 
nonsymmetric  terminal  sets  in  (2)  is  no  increase  in  generality. 

(3)  Let  Pf  be  a  measure  such  that  P9{X j  €  A)  —  Pq{X\-9  €  A),  and  P  —  Pq. 
Assume  that  X\  is  symmetric  under  P.  Let  P*  be  a  measure  under  which  the 
Wiener  process  has  drift  {.  Then 

(^  <  <}  -  P((IMio)l>  ,(‘o» 

+  ^{(*0  <  </»  <  '■)  +°(mT/s) 

Remarks. 

(1)  The  results  essentially  say  that  the  corrected  probabilities  for  the  random 
walks  are  the  probabilities  that  the  Wiener  process  crosses  the  boundary  pushed 


out  by  an  amount  equal  to  the  expected  excess  over  the  boundary.  This  has 
intuitive  appeal,  and  it  agrees  with  the  linear  cases  computed  by  Siegmund  and 
Yuh  [1982]. 


(2)  There  is  no  increase  in  generality  in  allowing  starting  times  other  than  0 
for  the  one-sided  boundary  problem.  This  can  be  done  by  conditioning  on  the 
position  of  the  random  walk  at  the  appropriate  time.  The  two-sided  problem 
requires  symmetric  boundaries  and  so  conditional  on  the  position  of  the  random 
walk  at  a  fixed  time  the  boundaries  appear  unsymmetric.  Evidently  the  symme¬ 
try  of  the  random  walk  makes  this  effect  cancel  out  on  the  kverage. 

(3)  There  are  many  variants  of  these  probabilities  that  it  would  be  nice  to 
do,  but  about  which  this  method  says  nothing.  Foremost  is  the  case  of  non¬ 
zero  skewness.  In  that  case,  this  method  will  derive  integral  equations  for  the 
correction  term  as  below,  but  no  solutions  can  be  exhibited,  except  in  the  case  of 
a  constant  boundary,  where  the  results  agree  with  Siegmund  and  Yuh's.  In  the 
two-sided  case,  it  would  be  nice  to  be  able  to  calculate 

as  well  as  to  be  able  to  handle  nonsymmetric  boundaries.  The  author  opines  that 
these  two  problems  are  linked,  and  that  there  is  an  “interaction”  term  between 
the  upper  and  lower  boundaries  which  makes  both  problems  difficult,  but  which 
cancels  out  if  sufficient  symmetry  b  present. 

(4)  In,  for  example,  a  normal  translation  family,  if  a  test  b  defined  by  a 

one-sided  boundary,  crossing  probabilities  for  contiguous  alternatives,  which  are 
those  whose  expected  value  b  follow  directly  from  these  results,  since  the 

drift  can  be  added  to  the  boundary.  However,  in  the  two-sided  case  the  result  b 
once  again  messed  up  by  the  introduction  of  asymmetry  into  the  problem.  The 
two-sided  case  b  the  more  interesting  one  for  applications. 

(5)  S.  V.  Nagaev  [1970]  establbhed  a  Berry-Essen  type  bound  of  tor 
thb  problem,  showing  that,  in  the  above  notation, 


where  L  is  an  absolute  constant,  by  using  basically  a  Fourier  inversion  method. 
That  analysis  is  very  difficult.  A  review  of  the  state  of  boundary  crossing  prob¬ 
lems  as  of  1972  can  be  found  in  Borovkov  [1972]. 


Heuristic  Proof  of  the  Result!  There  are  two  major  points  that  make  the 
proof  heuristic.  One  will  be  mentioned  when  it  arises  but  the  other  is  at  the 
heart  of  the  proof.  It  is  assumed  that 


P  {~  <  *}  ™  P{r  <  0 


where  C(t)  is  a  smooth  function.  An  integral  equation  is  then  derived  for  C(i), 
and  a  solution  to  the  integral  equation  is  exhibited. 

Consider  now  the  one-sided  boundary  problem,  so  that 


Tm  *  inf  {«*  :  Sn  >  ml!2f  (^)J. 


The  starting  point  for  this  analysis,  as  for  that  of  Siegmund  and  Yuh,  is  a 
simple  decomposition  of  all  paths  into  those  which  are  above  the  boundary  at 
time  I  and  those  which  are  not. 


-  p{^ji  >  *(«)}  +  "£  -  Sj  <  m'/2  (|(‘>-  »(£)) -R,Tm-  j] 

where  R  is  the  excess  over  the  boundary,  Sjm  ■■  m1/2  f(^)  +  R-  The  above  is 

p{^}i  > »<<>}  +  E>{^  -  si  <  ™1/2  (<*>  -»(£))  -  Rj}piT™  -  )}■ 
Now  Rj  e  a  random  variable  that  has  the  distribution  of  R  conditioned  on  {Tm  — 


m 


rr^wr^ 


j),  and  is  independent  of  Smt-Sj.  Now  by  an  Edgeworth  expansion  this  becomes 


Here  another  set  of  heuristic  assumptions  are  introduced.  First,  the  last  sum 
represents  a  quantity  of  the  form  ^jiE{f(Tm)R;  Tm  <  mt},  where  the  quantity 
f{Tm)  has  a  singularity  like  at  the  endpoint.  It  is  assumed  first 

of  all  that  the  distribution  oi  Tm  is  such  that  this  quantity  has  a  uniformly 
bounded  1  +  tth  moment  for  some  <  >  0.  Note  that  this  claim  is  true  for  the 
Wiener  process.  Siegmund  and  Yuh  also  ran  into  this  problem  ac  they  note 
in  [1082],  and  recently  Siegmund  has  circumvented  it  without  bounding  higher 
moments.  Nevertheless,  let  that  stand  as  one  assumption.  Second,  under  further 
assumptions  on  the  random  walk,  the  asymptotic  joint  distribution  of  R  and  Tm 
is  known  from  Section  1.  It  will  be  assumed  that  there  is  a  corresponding  local 
limit  theorem  which  says  that  ERj  —  ER  —  fi.  Combining  these  observations 
with  the  invariance  principle,  and  neglecting  terms  that  are  o  (^1/3),  the  following 
equation  for  the  hitting  probabilities  results: 


m 


1/2 


El  {TT7jl75 - 


One  further  simplification  is  possible.  One  sees,  upon  conditioning  on  TTf,  that 
the  last  expectation  Xdx  represents  P{W{t)  €  #(<)  +  dx,  rg  <  (}.  However 


jTJT—** 


P{W(t)  €  tit)  +  dx,  rg<t}  +  P{W(t)  €  tit)  +  dz,rg>t } 
P{Wit)€tit)  +  dx,  rg<t }. 


SO 


Therefore 


9i  0->(y/«n) 
y/t-H”* 


) 


\  m  mi  ml/2  (1/2 


(3.1) 


The  last  term  will  be  called  the  part  of  the  equation,  and  the  rest  the 

0(1)  part  of  the  equation.  This  is  to  be  regarded  as  an  equation  for  P  <  tj 
or  p|^  -  yj.  The  strategy  b  to  produce  a  function  that  satbfies  the  0(1)  part 
of  the  equation  except  for  terms  which  are  o  (^j),  and  to  represent  the  solution 
as  a  perturbation  of  thb  function.  The  candidate  b  f(m)  “  P 
To  see  that  thb  works,  consider  those  paths  for  which  rg  6 

increment  W  W(rg )  b  small  (E^W (4)  -  Mty))  *“  aQd  symmetrically 

dbtributed.  Therefore,  the  dbtribution  H,  defined  by 


rP€(*rir»  4)}- 

*ijr>  »£)•  The 


also  satbfies  EH 2  —  O(^),  and  EH  —  O(^).  Most  importantly,  conditioning  on 
rg  shows  that  f{m)  exactly  satbfies 


where  P(t)  —  /(^)-  A  simple  Taylor  expansion  now  shows  that  F(t)  satbfies 

the  0(1)  part  of  the  equation  except  for  terms  that  are 

Now  suppose  that  the  solution  to  (3.1)  can  be  written  as  F{t)  +  where 
Pj(l)  —  £mt  J\  Necessarily,  Fj(0)  —  0.  Then  an  equation  for  Pj(t)  can  be 
derived  by  substituting  thb  in  (3.1)  and  using  the  fact  that  F(t)  satbfies  the  0(1) 


- 


part  of  the  equation.  Explicitly 


*('(&  *-&*(£)) 

If  Fj  is  smooth,  /i(^)  -  and  neglecting  terms  that  are  0(1)  the  equation 
for  Fi(t)  is 


f(t)  -  g(j/m)\ 
y/t-  j/m  ) 

f  MVtW) 

~m>/2  ,1/2 


'(d^) 


As  an  aside,  note  that  integration  by  parts  converts  (3.2)  into  a  Fredholm 
equation  of  the  second  kind.  The  associated  homogeneous  equation  is 


(which  is  homogeneous  because  Fj(0)  »  0).  In  order  for  the  nonhomogeneous 
equation  to  have  a  unique  solution  it  is  necessary  and  sufficient  that  the  homo¬ 
geneous  equation  should  have  no  nontrivial  solution.  It  is  likely  that  this  is  the 
case,  however  the  solution  that  will  be  produced  has  sufficient  intuitive  appeal 
to  make  questions  about  its  uniqueness  almost  moot. 

Our  candidate  for  a  solution  is 


Substituting  this  in  (3.2)  shows  that  it  must  satisfy 

<  •>  u-.- 


There  should  be  no  problem  interchanging  orders  of  differentiation,  or  in  taking 
outside  the  integral,  and  so  doing  yields  for  the  right  hand  side  above 


(3.3) 


The  integral  has  a  clear  probabilistic  interpretation.  Consider  the  Wiener  process 
and  condition  on  7r then  clearly  the  integral  above  represents 


P{r0+h  <  *>  W(*)  <  t(t)  +  h)  -  P{r,+h  <t)-  >  rft)  +  h). 

And  so  the  expression  in  (3.3)  is 

4 '’('*+*  <  <)  l*_o 


~kk  akP' s+k  <  'a**0’ 


9h  t-p+n  -  *i  dh' 

'I 

~tH2  dh‘ 
which  is  to  say  that  &fiP{rg+h  <  0  U— 0  ^  a  solution. 

What  this  means  is  that 

p{7~  < '}  -  P{T°  <  ■ *MPlr***  <  •»  i‘-»  +-(m^) 

as  the  expression  in  the  first  line  is  the  Taylor  series  expansion  of  that  in  the 
second. 

The  second  statement  of  1  follows  directly  from  the  first  part: 


”  P  { <  *-  "W  <*('>-*  + 

-  <  *>  "W 6  (««)  -  *-  - 1 + + » (4/s) 

-  w(d  < ,«)  -  z) +.(-173). 


Now  consider  the  case  of  a  two-sided  symmetric  boundary  and  symmetric 
random  walk.  Here  0  <  fj  <  t  <  1,  Tm  —  inf(n  >  mti  :  |5n|>  m1/2  j(t)},  rg  — 


inf{<  >  «i  :  |W(0|>  m1/2  g(t)}.  The  analysis  that  leads  to  the  integral  equation 
is  very  similar  to  that  of  the  one-sided  case,  and  the  reason  for  the  restriction 
to  symmetric  situations  does  not  show  up  in  the  derivation,  except  insofar  as  it 
leads  to  a  different  integral  equation  to  which  the  author  can  present  no  solution, 
and  therefore  it  will  be  omitted.  The  equation  that  P  <  /j  should  satisfy  is 

P  {^T  <  *}  “  2(1 "  +  2(1 "  *(*(<2))) 

,  e  J  mo + MSI'  -  -  mo  -  «m /«  -o1'2) . \ 

+^mE( - - • '« < r  <  •) 

As  before,  some  analysis  is  possible  on  the  last  expression 

—  2P{W(t)  6  t(t)  +  dx,ti<r<  t,  W(*)  <  0} 

where  the  factor  of  2  comes  in  because  each  boundary  contributes  symmetrically, 
and 


Uitit )  -  f(r)/(t  -  r)1/2) 


“V  (t^r) V2  **  <  r  <  T* 

—  2P{W{t)  €  f(t)  +  dx,  tx  <  r  <  t,  W(r)  <  0}. 

Therefore,  the  sum  of  these  two  terms  is  2 P{W(t)  e  f(t)  +  dx,  t\  <  r  <  t),  and  so 
the  difference  can  be  represented  as 

4  P{W(t)  €  f(t)  +  dx,  W(r)  <  0,  «i  <  r  <  t)  -  2  P{W{t)  €  g{t)  +  dx,tl<r<t}. 
Call  this  H{t)dx. 

As  above,  the  Wiener  process  solution,  F0(/),  can  be  shown  to  satisfy  the  0(1) 
part  of  the  equation  up  to  terms  that  are  and  so,  if  a  solution  of  the 

form  Poit)  +  ^/2Pi(t)  is  posited,  it  is  easily  seen  that  fj(«)  should  satisfy 


The  proposed  solution  is  <  t,  |W(*1)j<  g{tx)}.  Substituting 

this  in  the  integral  gives 


which  is  the  ^  of 


P{rg+heds,  |W(t!)|<  ,(*!)}, 


P{<1  <  ',+A  <  «»  ^(rp+A)  >  °»  -#(«)  +  *<  MO  <  f(0  +  A,  !M*l)l<  *(*i)} 

+  P{*1  <  rg+h  <<»  M*^)  <  0, 

-  „(t)  -  A  <  <  y(t)  -  A,  |m«l)l<  *«1» 

-  ^{*1  <  rg+h  <  t>  Wi'g+h  >  0,  -1(0-  A  <  MO  <  f(0  +  A,  IM«l)l<  f(*l)} 

+p{rg+h  <  t,  w(rg+h)  <  o,  —g{t)—h  <  mo  <  ff(o+A,  im<i)i<  tih)) 

~P{h<rg+h<t,W(rg+h)>  0, 

-  g(t)  -  A  <  mO  <  -g(t)  +  A,  |m«l)l< 

P{tl  ^  rg+h  ^  ^(r^+A)  ^  ®» 

„(«)  -  a  <  mo  <  M + a,  im*i)i<  «(*i)} 

-  Pih  <  Tg+h  <  t,  |m*l)l<  *(0)} 

-p{imoi>  «(o+A,  <i  <  rg+h,  im«i)i<  ^i)> 

2P{t\  <  Tg+h  ^  t,  <  0, 

-  g(t)  -  a  <  mo  <  -*(0 + h,  im<i)i<  a(o)} 

-m  <  r9+h  <  t ,  im«i)i<  >«i)}-^{|moi>  j(o+a,  *i  <  t0) 

-  2P{tx  <  rg+h  <  t,  W(rg+h)  <  0, 

-  g{t)  -  A  <  mo  <  tit)  +  h,  |m*l)l<  IT(*i)}- 

The  ^  of  this  expression  is  easily  seen  to  be 

^P{t i  <  r*+fc  <  t,  Im«l)l<  U-0  +2P<mO€  lit)  +  dx,  tj  <  rp}/<fc 

-  4P{tj  <  fy  <  t,  W(Tg)  <  0,  mo  6  -p(0  +  dx}/dx, 
and  comparing  this  with  the  expression  for  H(t)  given  above  shows  that 

^P{r(+4  <  (,  WOK  ,<()> 

is  a  solution.  This  means  that 

p{lr < '}  ”  m‘/2 

+  P(«1  <  VfO/mW  <  '•  WlH<  »(<l))+»(^7j)- 
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The  second  statement  of  2  follows  as  did  the  second  statement  of  1. 

The  proof  of  the  third  result  is  a  little  bit  different.  Once  again  it  can  be 
argued  that  if 


then  satisfies  an  integral  equation.  The  exact  form  of  this  equation  b  not 
of  any  interest  except  that  it  must  be  noted  that  it  does  not  depend  in  any  way 
(other  than  symmetry)  on  the  particular  random  walk.  Now,  the  equation  itself 
constitutes  an  invariance  principle  for  the  correction  term,  and  if  the  correction 
term  can  be  found  for  any  member  of  the  family  by  any  method,  it  b  necessarily 
the  correction  term  for  other  symmetric  random  walks.  A  likelihood-ratio  method 
will  now  be  used  to  find  the  correction  term  for  normal  random  walk. 

Suppose,  therefore,  that  the  random  walk  is  Gaussian.  We  will  use  below,  in 
order  of  appearance,  Wald's  likelihood  ratio  identity  (Siegmund  [6],  propositions 
2.24  and  3.2)  on  the  random  walk,  Theorem  1.1,  result  2  of  this  section,  and  the 
likelihood  ratio  identity  on  the  resulting  diffusion  limits. 

'«/*/•  {£<•}-*  T”ln-  ~  <  <} 

-  e»  L«<w -»>+<>-<*/*  w«  £<(j 

-4"M!/2Wm(1 + £  <  •) 

-  *  W2m  (•  • ♦  &)  • £ <  •}  ■ ♦  ■ ’U) 

-  Eq  (Wm)K2/ 2  Tm/m.  Tm  <  A  +  #  f  M 

l  m  J  Vm1/2/ 


This  implies  the  stated  result. 
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